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Abstract 

We study dynamics of two coupled periodically driven oscillators. An 
important example of such a system is a dynamic vibration absorber which 
consists of a small mass attached to the primary vibrating system of a large 
mass. 

Periodic solutions of the approximate effective equation (derived in our 
earlier papers) are determined within the Krylov-Bogoliubov-Mitropolsky 
approach to compute the amplitude profiles A (SI). In the present paper 
we investigate metamorphoses of the function A (SI) induced by changes 
of the control parameters in the case of 1 : 3 resonances. 

1 Introduction 

In the present paper we analyse two coupled oscillators, one of which is driven 
by an external periodic force. An important example of such a system is a 
dynamic vibration absorber which consists of a mass m 2 , attached to the pri- 
mary vibrating system of mass mi [UH]. Equations describing dynamics of this 
system are of form: 



m\X\ — V\ (±i) — R\ (x\) + V 2 (&2 — x\) + R 2 {x% — x\) = / cos (ut) 

1Tl 2 X 2 - V 2 (#2 - Xl) - i? 2 (X2 — X\) = 



(1) 



where Vi, R\ and V%, R 2 represent (nonlinear) force of internal friction and 
(nonlinear) elastic restoring force for mass mi and mass m 2 , respectively. In 
the present paper we shall consider a simplified model: 

Ri Oi) = -ocixi, Vi (±i) = -vi±i. (2) 
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Dynamics of coupled periodically driven oscillators is very complicated, see 
[21 HJ |6j [7] and references therein. We simplified the problem described by 
equations (TTJ), deriving the exact fourth-order nonlinear equation for internal 
motion as well as approximate second-order effective equation in [5]. 

Moreover, applying the Krylov-Bogoliubov-Mitropolsky method to these 
equations we have computed the corresponding nonlinear resonances in the ef- 
fective equation (cf. [8] and [9] for the cases of 1 : 1 and 1 : 3 resonances, 
respectively). Dependence of the amplitude A of nonlinear resonances on the 
frequency u is significantly more complicated than in the case of Duffing oscil- 
lator and this leads to new nonlinear phenomena. In a recent paper we inves- 
tigated metamorphoses of the function A (ui) induced by changes of the control 
parameters in the case of 1 : 1 resonance [10] . In the present paper we continue 
this approach studying metamorphoses of A (w) for 1 : 3 resonance. 

In the next Section the exact 4th-order equation for the internal motion 
and approximate 2nd-order effective equations in non-dimensional form are de- 
scribed. In Section 3 amplitude profiles for 1 : 3 resonances are determined 
within the Krylov-Bogoliubov-Mitropolsky approach for the approximate 2nd- 
order effective equation (and for the Duffing equation which follows from the 
effective equation if some parameters are put equal zero). In Section 4 theory 
of algebraic curves is used to compute singular points of effective equation am- 
plitude profiles - metamorphoses of amplitude profiles occur in neighbourhoods 
of such points. In Section 5 examples of analytical and numerical computa- 
tions are presented for the Duffing equation. Our results are summarized and 
perspectives of further studies are described in the last Section. 

2 Exact equation for internal motion and its ap- 
proximations 

In new variables, x = x\, y = x% — x%, equations ([!]), ([2]) can be written as: 



where m = mi, m e = mj, v = v\, a = ai, V e = Va, R e = Rz- It is possible 
to simplify the problem eliminating the variable a; in (|3[) to obtain the exact 
fourth-order equation for the variable y only - describing relative motion of the 
mass m e [8]. 

In the present work we assume: 



m'x + vx + ax + V e (y) + R e (y) = / cos (tot) 
m e (x + y)- V e (y) - R e (y) = 





Re (y) = oi e y - j e y 3 , V e (y) = -v e y. 



(4) 



The exact equation for relative motion reads: 




(5) 
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where F = m e w 2 f , /i — mm e /M and A = m e /M is a nondimensional parameter 
i- 

Eqn. ((SJ can be written in the following nondimensional form [10] : 
t (0 + fc£ - z + z 3 ) + « (H & + a) = G^tf cos (fir) , 

where nondimensional time r and nondimensional displacement z of the mass 
m e are defined as: 

T = *D,z = yy /%, = (7) 

while nondimensional constants are given by: 

h = **r, H = t£=, n = ^. G = —xf^f, k = s±, a= -2^. (8) 

We shall consider hierarchy of approximate equations arising from ([B]) [10 . 
For small enough values of the parameters k, H, a we can reject the second term 
on the left in ([5]) obtaining the approximate equation which can be integrated 
partly to yield the effective equation: 



^+h^-z + z 3 = - 7 -, 2! cos(Or + 5), U^Grfa) (9) 

dr 2 dr y/(n 2 -a) 2 +H 2 n 2 v ; ' V K+1 J w 

where transient states have been omitted and tani5 = ?jrz- ■ And, finally, for 
H = 0, a = we get the Duffing equation: 

44 + k-T" - z + z 3 = -7 cos (fir + 5) . (10) 



3 Perturbation analysis of the 1 : 3 resonance 

The 1 : 3 resonance, a solution of the effective equation © of form z = 
A cos (3flT + ip), can be seen in the bifurcation diagram computed for the effec- 
tive equation - see Fig. 4.1 in [8], ui < 3.1. We apply the Krylov-Bogoliubov- 
Mitropolsky (KBM) perturbation approach [TTJ [12] to Eqn. (|9]) , working in the 
spirit of [6], to determine the corresponding amplitude profile, i.e. dependence 
of the amplitude A on frequency £7. 

To study subharmonic resonance 1 : 3 we cast equation ^ into form: 

+ 9 2 z + s ((-ao -B 2 + a z 2 ) z + h ¥) = -7" 2 co S (Qr+ A -) (n) 

dr 2 U u o u ; u dr ; ^/($2 2 -a) 2 +// 2 Q 2 v 7 

with 

eO 2 , = 9 2 , ea = 1, eh = h, (12) 
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where we assumed that the external force is of order e° rather than e 1 (see [T3] 
for discussion). 

We substitute z (r) = u (r) + uq (r) into (jTTJ) to remove the external forcing 
term on the right-hand side. We thus get: 

+ 6 2 u + ^ + & 2 u + eg (u, uq) = ~ 7 " 2 cos( " r+3) 



g (it, uq) = h d(u + T uo) + (-ao - ©o) ( M + u o) + a a (u + u f . 

Now we put uq (r) = C cos (Or + <5) into ([TU)) . It follows that for C = 

, „.) 1 „■> two terms on the left-hand side, 4ttt + © 2 un, and the 

x /(n2-a) 2 +/i-2n2 e 2 -n^ > dr- 

external forcing term on the right-hand side of (|13j) cancel out to yield: 

j^ + Q 2 u + eg{u,u ) -0. (14) 

We shall now determine approximate form of O 2 following procedure de- 
scribed in [6\. Neglecting in (jTTJ) the damping term /lo^f and external forcing 
we get: 

j^-z + z 3 =0. (15) 

Substituting in (fT5|> z (r) = A cos (6r), applying identity cos 3 (Or) = § cos (Or)-| 
j cos (30r), and rejecting term proportional to cos (30r) we get finally the ap- 
proximate expression O 2 = jA 2 — 1. 

We have thus written the effective equation © in form (fl"4")) with g(u,Uo) 
defined in (fT3"|) and: 

-o(r) = Ccos(Or + 5), ^7= ^ 7^^ , O 2 = \A 2 - 1. (16) 

Since we are looking for 1 : 3 resonances we have to consider frequencies Q 
close to 30. We thus put in (TT41 O 2 = (j) + ecr with cr of order e°, obtaining 
finally 

+ (§) 2 W + e(™ + 5 ( M , Uo )) = O. (17) 
We assume the following form of the solution: 

u = A cos (§r + ip) + eui (A,(p,r) + ... . (18) 

Substituting (|18l) into (1171) . eliminating secular terms and demanding that 
4^ = 0, ^ = to find stationary states we get finally [5]: 

, f o r\ -M? (19) 

tan ^ - S) = 3(f.4 2 + §C 2 -±0 2 -l)' 

with C given by (1161) . If we put H = 0, a = then we get implicit equation for 
the amplitude profile for the Duffing equation: 



m 2 + + \C 2 - ±0 2 l) 2 = (| AC) 2 , C = (3 _ A2 _l 2 _ iy (20) 
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4 Metamorphoses of the amplitude profiles for 
the 1 : 3 resonance 



Equations (TlTJl) , (|2"U1) define the corresponding amplitude profiles implicitly. Such 
amplitude profiles can be classified as planar algebraic curves, see [T3] for a 
general theory. Let L (X, Y; A) =0 defines such a curve where A is a parameter. 
A singular point (Xq,Yq) of the algebraic curve obeys conditions: 

i(W ) = 0, ^^)=0, ^™ = 0. (2!) 

Assume that a solution (Xq,Yq) of Eqns. (f2~Tj) exists for A = Ao and there 
are no other solutions in some neighbourhood of Ao. Let A < Ao, then the 
curve L (X,Y; A) = for growing values of A changes its form at A = Ao and, 
again, for A > A . We shall refer to such changes as metamorphoses (cf. [10] 
for metamorphoses of amplitude profiles in the case of 1 : 1 resonance in the 
effective equation). 

In the case of the effective equation the amplitude profile of 1 : 3 resonance 
is given by Eqn. (|19|) or, in new variables X = ft 2 , Y = A 2 , by the equation 
L (X, Y; a, 7, h,H)=0 where 

tt2- x2 



L (X, Y; a, 7, h, H) = U 4 (±h 2 X + U 2 ) ((X - a)' + H 2 X , J 
+3 7 2 U 2 X 2 (±Y -%X-1) ((X - a) 2 + H 2 X) + \^X\ ( 22 ) 



U = §Y — X — 1. 

Equations for singular points of the amplitude profile for the 1 : 3 resonance 
of the effective equation are given by (|21[) . (|2"2"j) . To find solutions of these 
equations we solve the following cubic equation: 

c 3 U 3 + c 2 U 2 + all + c = 0, 
c 3 = 567£ 2 

c 2 = -I8B1 (28A 3 + (87g - 81) X 2 + (146r - 54 g ) X - 27r) , , 

ci = -AB 2 X (967X 3 + 639A 2 + (162 - 967r + 423 9 ) X + 207r + 81q) ^ > 
c = 32B 2 X (55A 2 + 9) 

B 1 = 3X 2 + qX -r, B 2 = X 2 + qX + r, q = H 2 -2a, r = a 2 , 
for arbitrary A > 0, a > 0, H > and compute variables p, s: 



_ 2 (9U+8X)((3X 3 +2X 2 q+{-9q+r)X-18r)U-16X 3 -16X 2 q-16Xr) 
P ~ ~ 9 X(X81X 2 +27qX-27r)U-46X i > + (-18~4(>q)X 2 + {~46r-18q)X-18r) ' 



_ 4 U 2 (X 2 +qX+r)(l28X 2 + 162pX+360XU+243U 2 ) 
S ~ ~ 9 X 2 (46X+81l/-18) ' 

Finally, we find Y, h and 7 from definitions: 



(24) 



U = \Y-X-\, p=\h 2 , s = §7 2 , (25) 
and physically acceptable solutions must fulfill conditions: Y > 0, h > 0, 7 > 0. 
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We can obtain the case of the Duffing equation putting a = H = in the 
above formulae. We thus get: 



L (X, Y; 7, h) = C/ 4 (\h>X + U 2 ) + 3 7 2 [/ 2 ( ±Y - \X - l) + f 7 4 , (26) 

and 

5f03C/ 3 + (-1512* + 4374) £/ 2 + (-3868X 2 - 2556X - 648) U , 

+1760X 3 + 288X = ^ J 



where X is arbitrary, [/ = jY — X — 1 and 

i 2 _ f) 24[/ j y+27C/ 2 -128X-144L r 2 _ 8 r r 2 72LfX+243t/ 2 - 128X 2 fno\ 

,l Z 46X-8K7+18 ' ' 27^ 46X-8K7+18 ' ^ °> 



5 Analytical and numerical computations: the 
Duffing equation 

We shall find a metamorphosis of the bifurcation diagram for the Duffing equa- 
tion (jiop . To this end we have to compute a singular point of the amplitude 
profile flU- Let ft* = 1.6 (X» = 2.56). We get from Eqns. ([27]), O for 
X = X* one physical solution: U* = -2.946 246 654, h* = 0.894 933 8113, 
7* = 2.235 385 759, Y"* = 0.818 337 794 (A, = 0.904620248). 

In Fig. 1 we plot amplitude profiles, i.e. variables A, fulfilling (|2U)) . for 
the critical value 7 = 7* and h — 0.4, 0.8, 0.85 and the critical value h = h*. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 
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Figure 1: Amplitude profiles A(Q), 7 = 7*, h = 0.4 (blue), 0.8 (green), 0.85 
(sienna), /i* (red). 
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Figure 2: Bifurcation diagram for the Duffing equation, h = 0.4, 7 = 7*. 

Resonance f : 3 is shown in Fig. 2 where bifurcation diagram for the Duffing 
equation (fT0|) in (z, f2) plane was computed for ft, = 0.4, 7 = 7* and f2 € 
[2.3, 3.0]. 

Since the KBM method is approximate metamorphosis in the real system 
may happen at a slightly different value of, say, parameter h. The numeri- 
cally exact critical value of this parameter was determined from the bifurcation 
diagram below where dependence of z on h is shown for 51 = 57*, 7 = 7*. 



z 



0.5 




Figure 3: Bifurcation diagram for the Duffing equation, Q = 57*, 7 = 7*. 

It follows that the 1 : 3 resonance ends abruptly for growing h at h = 0.905, 
i.e. slightly above the critical value h* . 
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In Figs. 4, 5 below bifurcations diagrams showing dependence of z on ft for 
7 = 7* and h = 0.90, 0.92 are shown. 





Figure 4: Bifurcation diagram for the Duffing equation, h = 0.9, 7 = 7*. 

We realize that the 1 : 3 resonance disappears for growing h in agree- 
ment with analytical computations (based however on the approximate KBM 
method). 




Figure 5: Bifurcation diagram for the Duffing equation, h = 0.92, 7 = 7*. 
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6 Summary and discussion 



In this work we have studied metamorphoses of amplitude profiles for the 1 : 3 
resonances of the effective equation, describing approximately dynamics of two 
coupled periodically driven oscillators. Our analysis has been analytical al- 
though based on the approximate KBM method. Theory of algebraic curves 
has been used to compute singular points on amplitude profiles of the effective 
as well as the Duffing equation. It follows from general theory that metamor- 
phoses of amplitude profiles occur in neighbourhoods of such points. The results 
obtained can be compared with our work on metamorphoses of 1 : 1 resonances 
in the effective equation |10j . 

In Section 4 we have computed analytically positions of singular points for 
the amplitude profiles A (Q) determined within the Krylov-Bogoliubov-Mitropolsky 
approach for the approximate 2nd-order effective equation ([9]). In Section 5 an- 
alytical and numerical results have been presented for the case of the Duffing 
equation arising as the subsequent approximation of the effective equation. We 
have also computed numerically bifurcation diagrams in the neighbourhoods 
of singular points and indeed dynamics of the Duffing equation (1101) changes 
according to metamorphoses of the corresponding amplitude profiles. More ex- 
actly, we have found only the case of isolated singular point and this corresponds 
to creation or destruction of 1 : 3 resonance. We are going to investigate much 
more complicated case of the effective equation in our next paper. 

References 

[1] J. P. Den Hartog, Mechanical Vibrations (4th edition), Dover Publications, 
New York 1985. 

[2] S. S. Oueini, A. H. Nayfeh and J.R. Pratt, Arch. Appl. Mech. 69, 585 
(1999). 

[3] W. Szemplinska-Stupnicka, The Behavior of Non-linear Vibrating Systems, 
Kluver Academic Publishers, Dordrecht, 1990. 

[4] J. Awrejcewicz, Bifurcation and Chaos in Coupled Oscillators, World Sci- 
entific, New Jersey 1991. 

[5] J. Kozlowski, U. Parlitz and W. Lauterborn, Phys. Rev. E 51, 1861 (1995). 

[6] K. Janicki, W. Szemplinska-Stupnicka, J. Sound. Vibr. 180, 253 (1995). 

[7] A. P. Kuznetsov, N. V. Stankevich and L. V Turukina, Physica D 238, 
1203 (2009). 

[8] A. Oknihski and J. Kyziol, Differential Equations and Nonlinear Mechanics 
2006, Article ID 56146 (2006). 

[9] J. Kyziol, PhD Thesis, Kielce University of Technology, 2007. 



9 



[10] J. Kyziol, and A. Okniriski, Acta Phys. Polon. B 42, 2063 (2011); 
larXiv:1012.2T40V l[nlin.CD]. 

[11] A. H. Nayfeh, Introduction to Perturbation Techniques, John Wiley & Sons, 
New York 1981. 

[12] J. Awrejcewicz and V.A. Krysko, Introduction to Asymptotic Methods. 
Chapman and Hall (CRC Press), New York 2006. 

[13] N. N. Bogoliubov, Y. A. Mitropolsky, Asymptotic Methods in the theory of 
Nonlinear Oscillators, Gordon and Breach, New York 1961. 

[14] C. T. C. Wall, Singular Points of Plane Curves, Cambridge University 
Press, New York 2004. 



10 



